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Abstract 

Layered and pseudocubic Ru-based perovskites have been the subject of consider- 
able recent attention, due to their unusual magnetic properties and the discovery of 
superconductivity in one member of the family, Sr2Ru04. From a magnetic point 
of view, interest derives from the stable ferromagnetism in SrRuOa, gradually dis- 
appearing to a non-magnetic phase upon substituting Sr with isovalent Ca, a very 
unusual kind of behavior for 3d perovskites. On the superconducting side, inter- 
est was stimulated by theoretical conjectures and experimental indications that 
Sr2Ru04 might be a p-wave superconductor. We report first-principles LSDA cal- 
culations for ferromagnetic SrRuOs, antiferromagnetic Sr2YRu06, non-magnetic 
CaRuOs, and superconducting Sr2Ru04. In all cases, magnetic properties are well 
reproduced by the calculations. Anomalous properties are explained in terms of 
simple TB models and Stoner theory. An important result is that O bears sizable 
magnetic moments and plays an important role in the formation of the magnetic 
states. Based on these calculations, we have built a model for the g-dependent Stoner 
interaction, which we consequently applied to Sr2Ru04 to estimate superconduct- 
ing and mass-renormalization electron-paramagnon coupling constants. We found 
that spin-fluctuation induced p-wave superconductivity is possible in Sr2Ru04. The 
estimated critical temperature, specific heat and susceptibility renormalizations are 
all in good agreement with experiment. 



The recent discovery of superconductivity in the layered ruthenate, Sr2Ru04 
[1] has generated new interest in Ru-based perovskites. At first glance this ma- 
terial seems analogous to the high-Tc cuprates. For instance, it has a similar 
crystal structure (it is isostructural with La2Cu04) and is apparently close to 
a magnetic instability (Sra-Cai-j^RuOa and Sr2RuY06 are ferro- and antiferro- 
magnetic, respectively). On the other hand, the more we learn about ruthen- 
ates, the less similar to cuprates they seem. While initial interest was largely 
related to the similarity to the high-Tc materials, now it is more that ruthen- 
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ates are deemed interesting per se, and, at least in their magnetic properties 
they are more variegated and probably more interesting than cupratcs. Fm- 
thermore, it appears that superconductivity in Sr2Ru04 can hardly be under- 
stood without a good understanding of magnetism in ruthenium perovskites 
in general. Thus this paper naturally breaks into two parts. First, we discuss 
magnetism in ruthenates, specifically antiferromagnetic Sr2RuY06, ferromag- 
netic SrRuOs, and paramagnetic CaRuOs. We will show that despite the wide 
range of magnetic properties, they all are governed by a simple Stoner-type 
mechanism, which manifests itself differently depending on crystal structure. 
We then shall show how closeness to a ferromagnetic instability can produce 
a triplet superconductivity in Sr2Ru04 and explain its normal-state transport 
properties. We shall also discuss what is currently maybe the most intriguing 
question in the theory of superconductivity in Sr2Ru04, namely why the ex- 
periment shows finite electronic density of states at zero energy (in NMR and 
specific heat experiments) at as low as O.STc- 



0. 1 Magnetism 

The great majority of magnetic transition metal oxides are based on the 3d- 
series. Density functional theory in its standard local spin density approxima- 
tion (LSDA) does not work very well for some of these materials; it often fails 
to yield the correct magnetic ground state, in many cases it underestimates 
the magnetic moments, in some others it does not reproduce correct insulating 
behavior. In such cases it is customary to speak about "strong correlation be- 
havior" . The LSDA is essentially a mean field theory where electron-electron 
interactions are treated in an averaged way, and the nature of a magnetic in- 
stability is related to the standard Stoner model, where the paramagnetic sus- 
ceptibility, renormalized in the RPA-like manner, may diverge at some wave 
vector. On the other hand, in the strong correlation picture the zero order 
approximation is the large-f/ Hubbard Hamiltonian with an inherent antifer- 
romagnetic instability to it via the superexchange mechanism. The first thing 
to decide is which of the two basic approaches serves better as the starting 
approximation. 

An important mechanism for magnetic instabilities in a one-electron frame- 
work is the "Stoner model". This is a purely itinerant magnetism approach. 
In the LSDA the total energy is written as E = Tg + Eh + Ee_i + E^c, where 
Tg is the single-particle kinetic energy. Eh, Ef.-i, and E^c are the Hartree, the 
electron-ion, and the exchange- correlation energies, respectively. A ferromag- 
netic instability is, in this model, an instability with respect to a perturbation 
consisting of splitting the band by an exchange field, readjusting the Fermi 
level, and recalculating of E^c taking into account the created magnetic po- 
larization. It is easy to see that the energy between ferromagnetic and the 
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paramagnetic states in the lowest order in magnetization M is 

AE ^ —. (1) 

AN{0) A Sw? ^ ' 

The last variation, I — PExd^fri} ^ is called the Stoner parameter. It defines 
the renormalization of the paramagnetic susceptibility due to spin fluctua- 
tions, X = Xo/(l ~ -^Xo)- Note that when the exchange splitting is imposed 
upon a compound with more than one component, the total magnetization 
is expressed as M — J^iMi, where Mj is the magnetic moment of the i- 
th component and is proportional to its partial DOS at the Fermi level, 
Mi/M = Ni{0)/N{0). This lets one relate the average Stoner factor for a 
compound, I, with the Stoner factors of the constituent atoms: AAE^c = 
-HiMfli = -MY.i{Ni/Nfli, hence / = - Y.i{Ni/Nfli. Of course, actual 
LSDA calculations take into account distortions of the bands as a function of 
magnetization, as well as the higher order in M terms, neglected in the Stoner 
model. 

Looking at such ruthenates as Sr2RuY06, SrRuOs, and CaRuOs from the 
Stoner point of view, one observes that oxygen p-character is present at the 
Fermi level to a substantially greater extent than in the cuprates or most 3d ox- 
ides. Calculating the average / for these compounds one finds that the oxygen 
contribution, [Ai'o(0)/A^(0)]^/o is substantial; if it is neglected, the Stoner cri- 
terion IN{0) > 1 is not satisfied for any of them. If it is included, Sr2RuY06 
and SrRuOs appear to be unstable against ferromagnetic transitions, while 
CaRuOa, because of a slightly different DOS, is barely stable. Detailed anal- 
ysis of the magnetism in these compounds has been published elsewhere [2]. 
The key ingredient is the strong Ru-0 hybridization, which puts O character 
at Ep and assures the validity of the Stoner model. 

One can generalize Stoner approach to antiferromagnetic instabilities. The 
main difference from the ferromagnetic case is that the DOS in the Stoner 
formula has to be replaced by the one-electron susceptibility, N{0) = x(0) —>■ 
x{Q), where Q is the antiferromagnetic vector. What turns out to be impor- 
tant is that if the AFM ordering in question is such that some atoms do not 
bear a magnetic moment by symmetry, they should be excluded from the cal- 
culation of the average Stoner factor. This is the case in SrRuOs, and CaRuOs 
where oxygen, bridging two nearest neighbor Ru, cannot acquire a magnetic 
moment if the two Ru atoms are aligned antiferromagnetically. Correspond- 
ingly, the average / for antiferromagnetically ordered (Sr,Ca)Ru03 would be 
considerably smaller than for ferromagnetic analogues. In Sr2RuY06 there are 
no bridging oxygens and the ground state is antiferromagnetic, with the oxy- 
gens bearing a large fraction of the total magnetization. This is reproduced 
by detailed self consistent LSDA calculations. 
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0.2 Superconductivity 



LSDA calculations for Ru-based perovskites generally either predict a mag- 
netic ground state or a paramagnetic state very close to an instability. The 
quasi-2D Sr2Ru04 is not an exception — LSDA calculations give a Stoner 
renormalization (1 — NI)^^ = 9 (experiment gives similar numbers). Thus, 
one expects strong spin fluctuations to be present in this compound. The sit- 
uation is similar to Pd metal, where NI is also close to 1. It is very hard to 
expect that a conventional superconducting state would survive in the presence 
of such spin fluctuations. In fact, Pd has a sizable electron-phonon interaction 
and would have been a superconductor apart from spin fluctuations, and in fact 
becomes such in amorphous state where spin fluctuations are suppressed [3]. 
On the other hand, it is known (see, e.g., Ref.[4]) that spin fluctuations provide 
effective repulsion for the singlet (s, d) pairing, but attraction for triplet (p) 
pairing. Thus it is tempting to ascribe superconductivity in Sr2Ru04 to the 
spin-fluctuation induced p-wave pairing[5,6]. LSDA calculations can be used 
as a tool to get a feeling about the size of the attraction provided by exchange 
of spin fluctuations and whether it is sufficient to explain the superconducting 
and normal state properties of this material. 

The valence bands of Sr2Ru04 are formed by the three t2g Ru orbitals, xy, yz, 
and zx. These are hybridized with the in-plane oxygen and, to a considerably 
lesser extent, with the apical oxygen[7,8] p-states. The bare oxygen p levels 
are well (~ 2 eV) removed from Ep, so the effect of the O p orbital is chiefly 
renormalization of the Ru t2g levels, and assisting in the d — d hopping. With 
nearest neighbors only, this gives one nearly circular cylindrical electronic sheet 
(7) of the Fermi surface (FS) and four crossing planes (quasi-lD FS). The weak 
xz — yz hybridization reconnects these planes to form two tetragonal prisms, a 
hole one (a) and an electron one (/3). De Haas- van Alphen experiments conflrm 
this fermiology[9]. In fact, the LDA a, f5, and 7 areas deviate from the dHvA 
experiment by only -2%, -3% and 5% of the Brillouin zone area, respectively, 
and an exact match can be achieved by very slight shifts of the bands a, f3, and 
7 by 5, -4, and -3 mRy, respectively. Such agreement is generally considered 
very good even in simple metals, and the small mismatch (which does not 
change the FS topology) is may be due to some underestimation in LDA 
calculations of the tiny xz—yz hybridization. Both calculation and experiment 
give nearly two dimensional Fermi surface: the relative c-axis variation of the 
extremal cross-section areas of the sheets /? and 7 is 6% and 1.5%, respectively 
(for these two sheets the extremal cross-sections are in the planes kz — 
and kz = tt/c). For the sheet a the relative change is 2% (for this sheet the 
extremal cross-sections are in the planes kz = and kz = 7r/2c). Experiment 
gives the numbers about twice smaller for all three sheets[10]; the difference 
is larger than the computational error, and presumably has its origin in the 
effects beyond the density functional theory. We repeated the calculations 
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using two non-LDA techniques, generalized gradient approximation[ll] and 
weighted density approximation [12], but the numbers hardly changed. In the 
following all calculational results are from the LDA linearized augmented plane 
wave calculations [7]. 

We assume that the exchange of the spin fluctuations is responsible for super- 
conductivity (and for the mass renormalization, to be discussed later). Such 
an interaction in metals was studied with respect to possible superconductiv- 
ity in Pd in the late 1970-ties (see, e.g., [13, 14]), and later in connection with 
heavy fermions. Assuming the Migdal theorem (a common approximation, al- 
though not well justified for spin fluctuations), the parallel-spin interaction, 
relevant for triplet pairing is given by the sum of the bubble diagrams with 
odd numbers of loops. 

Via -I, I,')- -^'(g)Xo(g) 



Here xo is the one-electron susceptibility, given as 

/kg ~ /k+q„fl 
ka/3 ^ka ~ ^k+q,/3 



Xo(q) = J2 — — exp(iqr) |k + q,pf , (3) 



with the usual notations. We used the approximation[15] Xo(q) = Xo(0) = 
N{0); this is a good approximation for an isotropic two-dimensional Fermi 
liquid[16]; we are currently investigating the quality of this approximation 
for Sr2Ru04, which is a two-dimensional, but not isotropic, Fermi liquid (so 
some modification of x(g) due to Fermi surface nesting may be expected). 
In any case, the g-dependence of I{q) is to be taken into account. As dis- 
cussed in the previous section, for the antiferromagnetic arrangement Iafm = 
I {TV /a, TT/a) = Iru {NRu/Nf, while Ifm = /(O) = Iru {NR^/Nf+2Io {No/2Nf . 
Atomic Stoner factors for Ru and O ions are calculated in a standard way and 
arc Iru ~ 0.7 eV, Iq ~ 1-6 eV. We found Iafm to be smaller than IpM by 
14% (oxygen contribution A/ = 0.06 eV). A g-dependence that reflects this 
effect is I{q) = //(I + b\^), where 6^ = 0.5(a/7r)2A//(/ - A/) « 0.08{a/TTf. 

Using these numbers, we calculate the effective coupling constant in p-channel. 
Following the suggestion of Agterberg et a/[17], we calculate the coupling con- 
stants separately for the three bands in question: xy (7), yz {(), and zx (^). 
The corresponding formula is 

A?, = {N,Nj/N){V(k - k')(vL • K)/(vlvi,)),j, (4) 



where i and j label the three bands, and v is the Fermi velocity. By symmetry. 
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the coupling matrix is 



( \p \P \P \ 
^S7 ^S? 







(5) 



and we calculate A?J^ = 0.16, A^^ = 0.075, and A^^ = 0.025. The critical tem- 
perature is defined by the maximum eigenvalue of the matrix (iV/iVj)Afj [18]. 
The corresponding eigenvector defines the relative magnitude of the order pa- 
rameter in bands 7 and () near Tc. We find the maximum eigenvalue of the 
corresponding coupling matrix is Xp =0.43, and the corresponding supercon- 
ducting state is O.857 -|- 0.38.^ -|- 0.38^. It is worth noting that using notations 
of Ref.[17], and taking into account the partial DOS : : ^ 0.44 : 
0.28 : 0.28, the matrix (5) can be translated to the interaction matrix U as 



U 



/ \ 

Vj<y^ Vj^Oi '^^13 

U'yci Uaa 



(6) 



where u^^ : Uja '■ '■ Uaa '■ Ua/s '■ UfSfs = 0.96 : 0.08 : 0.16 : 0.25 : 0.51 : 1, to be 
compared with the value conjectured in Ref.[17], 0.09 : 0.09 : 0.09 : 1 : 1 : 1. 
Their hypothesis about the smallness of the nondiagonal elements Uja and it-y^j 
is confirmed by the calculations, but the assumption about the smallness of u^^ 
is not. In any event, the calculated value of Xp =0.43 is sizable, and STifRcicnt to 
explain the observed superconductivity. We would like to emphasize the role of 
oxygen in this scenario: if not for the oxygen Stoner factor, the g-dependence 
of the effective interaction V{q) would be so small that the coupling in Eq. 
(4) would average near zero. 



0.3 Renormalization 



The mass renormalization is not as easy to define. Besides the parallel-spin 
interaction (2), there is the antiparallel-spin interaction, given in the same ap- 
proximation by the sum of the chain diagrams with even numbers of loops, plus 
ladder diagrams [13,19]. In the case of a contact interaction, the total inter- 
action is three times stronger than the interaction in the parallel-spin channel 
only. It was pointed out [14], though, that there is no good physical reason to 
single out any particular class of diagrams. It was found that including all three 
classes above leads to systematic overestimation of mass renormalizations by a 
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factor of 2 to 3 [13,20]. The present case is further complicated because unhke 
the electron-phonon interaction, the electron-electron (and. correspondingly, 
the electron-paramagnon) interaction is already included in some average way 
in the LSDA band structure. Thus, the electron-paramagnon mass renormal- 
ization is to some extent included in the LDA mass as well. 



Despite all these difficulties, one can get an idea about the size of the electron- 
paramagnon mass renormalization by making calculations with the parallel- 
spin interaction (2) only. The mass renormalization then is computed in the 
same way as the electron-phonon renormalization, i.e, by taking the average 
of V{q) of Eq.(2) over the FS. One has to remember, though, that there are 
other effects beyond the LDA, apart from the one that we calculate, which 
may further increase the observable mass. 

The coupling matrix which defines mass renormalization is written as A^j — 
(NiNj/N) < V(k — k') >ij, and the mass renormalization in band i is de- 
fined as A| = i^i^'^ J2j Aij. The average mass renormalization is A* = J^ijA^j- 
We calculate A;^ = 0.35, A|^ = 0.32, A^^ = 0.16, A|^ = 0.03. This gives 

a; = (a;^ + 2A;^)/i/^ = 1.5, A| = (A|^ + a;^ + ai^)///^ = 1.8, a^ = 1.7, to be 

compared with experimental dHvA values of 3, 2.3, and 3, respectively. The 
difference may be due to an electron-phonon coupling of the order of 1 and/or 
antiparallel spin fiuctuations, neglected in our calculations, as well as to the 
omission of the non-Migdal diagrams. In view of the underlying approxima- 
tions, the agreement is fairly good. 

One of the key problems, as discussed in Refs. [21,17], is the residual electronic 
specific heat [22], which remains at about 50% of its normal value well into 
the superconducting regime. There are superconducting solutions ("nonuni- 
tary states") for triplet pairing that are gapless, that is, have finite density 
of states at zero energy and zero temperature. However, the pairing energy 
for such states is lower than for the gapped states considered above. This 
led Agterberg et al [17] to postulate a pairing matrix that yields a vanish- 
ing gap for the 7 band. This, however, does not square with the quantitative 
estimate presented here. An earlier assumption[21,6] was that the excess pair- 
ing energy that forbids nonunitary combination of the order parameters may 
be overcome by additional magnetic (Stoner) energy in a nonunitary state. 
The requirements are strong Stoner renormalization (supported by the cal- 
culations) and strong particle hole asymmetry[23]. However, a quantitative 
estimate according to Ref.[23] shows that the effect is by far too weak. The 



criterion is 



Ted log Af 
dEp 



I- IN ^ ~ 10 ^, while it should be of the order 1 



for the nonunitary state to exist. 

Another possibility is related to an observation made a decade ago in connec- 
tion with the high-Tc superconductivity[24]: A well-known fact is that virtual 
phonons, even in a strongly coupled system, have no pair-breaking effect, so 
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that the density of states remains zero below the gap at zero temperature in 
a clean superconductor. However, this is a consequence of an internal symme- 
try of the Eliashberg equations, namely that the coupling function Q;^F(a;), 
entering the equation on A, is the same as a^F{u), entering the equation 
on Z. In case of p-wave pairing, for instance, this is not true any more, and 
formally there is finite density of states inside the gap at any temperature. 
Unfortunately, direct calculations [25] show that this effect is quantitatively 
strong only if a noticeable part of a'^F^cu) exists at a; < A, which is not the 
case in Sr2Ru04. 

Maybe the simplest explanation of the "residual DOS mystery" is still the 
most plausible. Despite the large mean free path, which in the reported 1.35 K 
samples reaches 1500-2000 A[26], this superconductor is still in the dirty limit: 
the Abrikosov-Gor'kov pair- breaking parameter 7 = 1/2tA = 7r^o/2im.f.p. = 
0.7, using the value for the coherence length ^0 — 1000 A[27]. Nonmagnetic 
impurities in a unitary 2D p-wave superconductor act as magnetic impurities 
in an s-wave superconductor. The DOS is given by the standard expression 

where u{E) satisfies the equation 

E = u — ^xj \/l — 

The resulting DOS at T = Tc/3 is shown on Fig.l and is seen to be very large 
below the gap (and does not show any trace of piling of the DOS above the 
gap). 

0.4 Conclusions 

To summarize, we have presented first principles calculations indicating that 
interactions due to exchange of FM spin fiuctuations, as calculated from the 
LDA band structure, are sufficiently strong to explain both the mass renor- 
malization and superconducting critical temperature of Sr2Ru04. 

This work was supported by the ONR. Computations were performed at the 
DoD HPCMO NAVO and ASC facilities. 
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Fig. 1. Relative density of states at T = O.STc in Abrikosov-Gor'kov theory for 

pair-breaking parameters 7 =0, 0.07, and 0.7. We estimate that for 1.35 K super- 
conducting samples the pair breaking parameter 7 is at least 0.7. 
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